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Abstract 

The statistical mechanics of spin models, such as the Ising or Potts models, on generic random 
graphs can be formulated economically by considering the N ^ 1 limit of N x N Hermitian matrix 
models. In this paper we consider the A'^ ^ 1 limit in complex matrix models, which describes vertex 
models of different sorts living on random graphs. From the graph theoretic perspective one is using 
matrix model and field theory inspired methods to count various classes of directed graphs. 

We also make some remarks on vertex models on planar random graphs (the N —f oo limit) where 
the resulting matrix models are not generally soluble using currently known methods. Nonetheless, 
some particular cases may be mapped onto known models and hence solved. 



1 Introduction 



In recent papers we have discussed various analytical and numerical aspects of an economical way of 
describing spin models living on "thin" (i.e. generic, non-planar) random graphs 0, In particular, one 
finds mean-field like behaviour for the spin models due to the locally-tree-like structure of the graphs. 
This approach is based on the simple observation [|| that the thin graphs appear as the Feynman diagrams 
in the perturbative expansion of matrix models when the ^ 1 (scalar) limit is taken. An alternative 
iV — > oo limit is perhaps more familiar for matrix models in the context of two-dimensional gravity, where 
the resulting "fat" graphs are of interest because of their relation to surfaces and string worldsheets. This 
relation is lost in the thin graph case |^ but the models still merit attention, both from the point of view 
of investigating various decorated random graphs and as a way of accessing mean field behaviour without 
infinite range interactions or the boundary problems of genuine trees like the Bethe lattice. 

The thin graph models are amenable to a saddle point solution, using large order/instanton methods 
from path integrals Q for the ordinary integrals which appear in this case. Such saddle point methods 
in graph theory been independently derived from a probabilistic viewpoint by Whittle in The topic 
of the current paper has also been presaged by the same author in where the main concern was 
understanding the statistics of random directed graphs using complex integrals rather than the real 
integrals of ||]. We take a rather different tack here, where our aim is to investigate vertex models 
on random graphs in their own right, so the focus is on the matter living on the graphs rather than the 
graphs themselves. 

In diagrammatic terms the effect of going from a Hermitian to a complex matrix model is to place 
arrows on the edges of the graphs that appear in the perturbative expansion, giving a directed graph. 
We take (/>''' to be the "head" of the arrow by convention. This is true whatever the value of N so it still 
holds, in particular, for the iV — > 1 limit. We can pick out the 2nth order in the perturbative expansion 
by carrying out a contour integral in the coupling constant in an exactly analogous manner to the path 
integral calculations in [Q . The general form of the partition function we are interested in for a vertex 
model on an ensemble of thin random graphs with 2n vertices is thus 
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where (f> is complex, (j)^ is its conjugate, P{<j>\ (f>, A, {A}) is the potential that specifies the particular 
vertex model, A is the overall vertex coupling that will invariably be scaled out in the ensuing discussion 
and {A} is the set of remaining couplings that can be tuned to search for critical behaviour. The prefactor 
Nn counts the number of undecorated graphs in the class of graphs of interest and will generically grow 
factorially with the number of vertices. For (j)'^ (3-regular) random graphs, for instance, it is given by 



N.„ 



(6n- 1)!! 
(2n)!! ■ 



(2) 



It is well known that complex matrix models with potentials of the form P{{(j>''(f)), A, {A}) generate 
"chequered" surfaces in the N ^ oo limit, where all the loops have an even number of edges [Q. In this 
case the critical behaviour turns out to be identical to that of a Hermitian model, as restricting oneself to 
combinations of effectively gives a Hermitian model in terms of Something similar occurs in the 
thin — > 1 limit, as was noted already in When one considers models where the in-degree (number 
of arrows in) and out-degree (number of arrows out) at each vertex is the same, as would be the case for 
potentials of the form P(((/)t0), A, |A}) the saddle point equations are singular in the original variables 
due to the radial degree of freedom q Switching to polar co-ordinates eliminates the singularity and one 
recovers the critical behaviour of the equivalent real integrals. 

The natural symmetry to preserve in the vertex models is arrow reversal, or complex conjugation. 
In the integral for the partition function in equ.(^ this corresponds to demanding a real integrand. We 

^The surfaces that might be associated with the graphs are "infinite genus", with a maximal handle density. 
^An identical situation occurs in the path integral for coupled anharmonic oscillators, where the large orders saddle point 
solution has a divergence at the radially symmetric point [M. 
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now note that we are not obliged to restrict ourselves to potentials composed only of combinations of 
to maintain this symmetry, which allows us to look at a more general class of potentials and observe 
a much richer critical behaviour than that seen for P{{(j)''(j)), A, {A}). We shall see in the sequel that 
the resulting models play a similar role to the eight-vertex model on a square lattice in that they contain 
many other soluble models as special cases. 

In the next two sections of the paper we discuss models on thin (f>^ and (f>^ graphs separately. In all 
cases we shall convert our models back to a real notation by making the substitution cj) — x+iy, (f)^ — x~iy 
{x,y real), as well as pre-emptively scaling out the overall vertex coupling A in the actions we discuss. 
As it is the critical behaviour of the matter that is of interest, which is determined by the saddle point 
equations for the "fields" ((/)^, (f), we shall concentrate exclusively on them, omitting the trivial saddle 
point in A. We then move on to discuss vertex models on planar graphs (i.e. the N oo limit), where the 
resulting matrix models arc not soluble in general, highlighting some cases where solutions are available 
by mapping onto previously solved models. 



2 Vertex Models on (Thin) 03 Graphs 

On 0"^ graphs we can have vertices with 3 inward or 3 outward pointing arrows, 2 arrows in and 1 out, 
or 2 arrows out and one in as shown in Fig.l. These correspond to terms of the form (0^)'^, (j)^ , 
and (p^icj))^ respectively in P{(t>\ (j), A, {A}). We can construct four possible linear combinations of these 
that are invariant under the complex conjugation symmetry, namely 

[(0t)20 + ^2(^t)]^ Z[(0t)20_^2(^t)] (3) 

SO the most general invariant action S {— ^(j>''4' + P{4>\ <j>, A = 1, {A})) one can write down is 



or, grouping together the similar vertices 

S = i |<^t^ „ (^^(^t)3 _ (^ + M^3 _ + [(</>t)^^] - (7 - ^5) [^^(</.t)] I . (5) 

We thus see it is natural to have the , (j>^ and (fy^icj)^) vertices appearing with complex (but 

conjugate) weights. We can cast this action in real form by writing (j> = x + iy, (p^ — x — iy, with x, y 
real which gives 

S = lix' + y') - (^±M^3 _ (^±My3 _ _ ^)^y2 _ _ ^)^2y^ (6) 

This transformation is rather reminiscent of the line decomposition that is used for vertex models on 
regular lattices with a; being black and y white lines, but we do not have line conservation at each 
vertex for the jjiodel. The saddle point equations dS/dx, dS/dy = for this action are explicitly 
soluble for general values of the couplings, but the solutions are not terribly illuminating for the purposes 
of extracting different critical behaviours when written down baldly. 

We can get a little more insight by reducing various subsets of the coupling constant space to already 
known models. One obvious reduction is to set P = 5 — which removes the terms antisymmetric in 
(f)) in the action. This leaves one with 

S=\[x'+v')-''^^x''[^~a)xy' (7) 
which displays mean-field Ising critical behaviour on a suitable coupling constant locus. 
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This can be seen by considering the standard action for the Ising model on graphs 

S^\{x^+V^)-cxv-^-{x^ + y^) (8) 
(where c = exp(— 2/3)) whose Feynman diagram expansion represents Ising spins with Hamiltonian 

H = (5Y,a,a, (9) 

Hving on the graphs, where the spin interactions are along nearest neighbour < ij > edges. If we carry 
out the linear transformation a; — > (a; + y) , y ^ -^[x — y) and some judicious rescalings we obtain 

S^^{x'+y^)-^[x' + 3c*xy'], (10) 

where c* — (1 — c)/(l + c), which is identical (up to a rescaling) with equ.(|^), provided we identify c* 
with (7 — a) /{a + 87). 

The Ising model on (j)^ graphs described by the action in equ.(^) displays a mean field transition 
at c = 1/3, so the model in equ.(p^) would be expected to display a transition at the "dual" value of 
c* — 1/2. A direct solution of the saddle point equations for equ.(|lC|) 




^ (mghT*), (11) 

where the T* is the temperature associated with the dual coupling, c* — exp(— 2/3*), reveals that this is 
indeed the case. It is worth remarking that the action of equ.(0, |lO|) on planar graphs (i.e. in the N 00 
limit) corresponds to an Ising model with spins situated at the vertices of a random triangulation dual to 
the original 0'^ graph |l^. The notion of duality is not defined for the non-planar graphs of the ^ 1 
limit but one could think of action representing spins residing on the vertices of the triangles in the cactus 
(or Husimi tree) diagram in Fig. 2 interacting with their nearest neighbours. This cactus graph is the 
medial graph of the underlying (j)^ graph which is also indicated in the figure. Unlike standard cacti, 
the branches of triangles eventually close to form large loops, since the underlying (j>^ graphs are closed. 
For an alternative physical interpretation one could consider the (j>^ graphs themselves and think of the 
model as a gas of y loops in an x "sea", which is the 0(1) loop gas representation of the Ising model. 

Viewed as Ising spins on a cactus graph the action of equ. (0, |l^) restricts triangles to have all spins 
up {x^ vertices) or two spins down and one up (xy^ vertices). This restriction does not appear in the 
planar case as one can then think of x representing a "same-spin" edge and y a "different-spin" edge, 
which is sufficient to cover all configurations when the triangles are pasted together along their edges. 
The difference on the thin graphs is that the triangles are to be thought of as being joined at their 
vertices. Lifting the constraint is equivalent to allowing non-zero S and (3 in the original vertex model 
and we obtain actions that describe multi-spin interaction models on cacti. The spin Hamiltonian for 
such a model is given by 

A i 

where the three spin sum is over triangles A. With all four possible spin configurations on the triangles 
of the Husimi tree now allowed, namely (-1- + +), (— — — ), {+ + — ), (-1- — — ), the action can be 
written as 

^ ^ 2^ + 2^ - a; y - C/^y a; - —x - -y , (13) 

where ^ — exp(2/3/i), ( — exp(2/3J3). This is clearly equivalent to equ.(^) once one rescales x x/^/JI 
and makes the identifications 

P + 3S = 1 
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a + 37 = Cm^/^ 
7 - a = Cm^^^ 

<5-/3 = M (14) 

which gives two conditions on our original four vertex model couplings to match the number of physical 
parameters. Note that it is natural to include the external field h from the start in such models as the 
phase transition generically appears in non-zero field for such multi-site interaction models Jill . 
The Hamiltonian of equ.(|2|) may be extended to include nearest neighbour interactions as well 

aiajCTk + J2'^cria-j + h'^ai). (15) 

A ij i 

and a similar reasoning to the pure 3-spin interaction case above leads to the action 

S = V {^^x^ + ^?x'y - Cuv'x - - ^y^) , (16) 

where the additional parameter v = exp(2/3J2)- Once again a rescaling x x/ ^JJi and absorbing the 
overall factor of v in the vertex coupling gives an action of the form equ.(^). We can also see that sending 
J3 ^ oo (C ^ 00) picks out the two sorts of vertex/triangle allowed in the Ising model discussed above. 
From the loop gas point of view both equ.( |l3| , p^ represent a gas of y loops and lines in an x sea, which 
can now interact via y^ terms and break via x^y terms. Alternatively the roles of y (loops/lines) and x 
(sea) may be exchanged as their vertices appear in a similar form in the action. 

It is also possible to bring Potts models under the umbrella of the vertex model framework in one of 
two possible ways. The original q-statc Potts model action contains q (real) variables (pi 

S^\jZ^l-cY,<P^^,~\jZ't'l (17) 

where c is now l/(exp(2/3) + q — 2), and represents the fohowing HamUtonian on graphs 

i7 = /3 ^(<5,.,,^. -1), (18) 

where the spins ai take on q values and 5 is a Kronecker delta. It is possible to write an "effective" action 
coming from imposing the observed symmetry breaking pattern 4>i...q-i = x,(f>q — y from the saddle 
point equation solutions on the original q variable action Q. This effective action 

S=\{q- 1) [1 - c{q 2)] x^ - i(g - 1):^^ + ly^ _ 1^3 _ _ ^^^^ (^g) 

faithfully reproduces the correct saddle point solutions. One can carry now out some rescalings on 



equ.(19) along with the ubiquitous linear transformation a; — s- (x — y)/\/2, y — > (x -|- y)/^/2 to obtain 



2' " ' 3 

where 



S = hx"^ +y'^)- \ [x^^ + SflAx^y + 3n^xy^ + Q^Ay^] (20) 



X = (9-l)(l-c((7-2)) , K^c{q-1) 
(q-D- y3/2 
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which is yet another variant of the action in equ.(^ for a particular choice of couphngs. It is noteworthy 
that the first order transition point for the g-state Potts model that was derived in [Q 

^ = ^2 ^22) 

is the point at which A = 0, reducing equ.(^o|) to an Ising-like action . 

A second possibility is to include two further quadratic terms in the allowed action 

^[02-(0t)2] ^ _2xy. (23) 

which are best thought of as addition bivalent vertices from the point of view of maintaining a vertex 
model interpretation. These allow two arrow heads or tails to meet. With these additional vertices the 
action becomes 

- ^^^{c^' f - ^^±^<A3 - (7 + ^5) [(^t)2^] _ _ ,5) f^2(^t)] I . (24) 

and choosing a = 7 ^ (g - l)/4, P = 5 = 1/4, ^ = -c(g - l)/2, ^ = [I + {q - 1)(1 - c(g - 2) )]/2, 
V = [\ — {q~ 1)(1 — c{q — 2) )]/4 gives an action which is recognisable as that for a q-state Potts model 
in equ.(p9|) when written in real form. 

To summarize the results of this section - the action of equ.(||) for a vertex model on thin (j)^ graphs 
has been shown to contain as special cases Ising and multi-spin interaction models, results which are 
easiest to see once the model has been rewritten in the real form of equ.(|^). Written is this form the 
vertex models can be re-interpreted as representing various systems, including spins on cacti or loop/line 
gases on the original ^'^ graphs. Potts models may also be incorporated into the general framework, either 
by reconstructing the effective Potts model action by rescaling and linear transformations from equ.(^, 
or by introducing further quadratic vertices as in equ.(^). 



3 Vertex Models on (Thin) 0^ Graphs 

We deal with the (/)^ vertex models rather more briefly, as the essential features are similar to the (j)^ 
models in the previous section. For (f)^ graphs one has the option of being able to preserve the in-degree 
and out-degree of each sort of vertex, but as we have noted this is equivalent to considering potentials 
composed of products of (f^ (f) and exploring a rather restricted class of critical behaviours. In a similar 
manner to the (jy^ graphs, we therefore consider all the conjugation symmetric terms 

[(</)t)4 + 04], ^[(04_(^t)4] 

[(0t)3^ + 03(^t)]^ Z[(0t)3^_^3(^t)]. (25) 

The most general action containing all these terms is 

S ^ l|<^V„|[(0t)4 + ^4]_^^[^4_(^t)4]_^(^t)2^2 

- a[(^t)3^ + ^3(^t)] _,^[(^t)3^_^3(^t)]} (26) 

or 

S=\ - (^±M^4 „ (^_M(^t)4 _ ^(^t)2^2 _ (^ ^ ^p)m^■^ _ (^ _ ^P)[<|>H^^)]] ■ (27) 
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As for the examples, this is best digested in bite-size chunks. Switching to real notation 

^ + ,2) _ + 27 + 4.) ^4 _ (c^ + 27 - 4.) (3. + 2^) 

2^ ' 4 4 2 

- (2p + /3)x3y-(2p-/3)a;y3 (28) 

we again find an action that has Ising and multi-spin interaction critical behaviour for various reduced 
subsets of its coupling constants. Setting the asymmetric terms p and [3 to zero again gives mean-field 
Ising critical behaviour, as can be seen in an exactly analogous fashion to the 0'^ graphs by carrying out 
the linear transformation x ^ -^(x -)~ y), x — > '^{^ ~ v) the action for an Ising model living on 
graphs 

to get (with suitable scalings) 



,4 



S=\{x^+y')~cxy~\{x^ + y^) (29) 



c 1/2, 2\ I4 3^22 1/ *\2 4 ior\\ 

^^2^^ +y -^cxy --{c) y . (30) 

The original (p^ action displays a mean-field Ising phase transition at c = 1/2, so one would expect the 
transformed action to have a transition at c* — 1/3, which a direct solution of the saddle point equations 
for equ.(|ol) 

a; = 1, y = {LowT*) 



confirms, just as for the (j)^ Ising model. 

Once again the Ising- like actions are seen to reside in the symmetric subspace of couplings. Allowing 
non-zero p and /? in equs.( [2^j2^ ) gives actions that may be interpreted as representing either multi-spin 
interactions on cacti composed of squares, such as that shown in Fig. 3, or loop/line gases on the original 
(j)'^ graphs. The spin model Hamiltonians on cacti graphs represented by such actions are of the form 

O-iCTjO-feO-; + J2 ^ CTiCTj + ft.^ (Tl), (32) 

[i,],k,l] ij i 

where the sum [i,j,k,l] is round the squares of the cactus, with similar relations between the spin 
couplings ( J4, J2, h, /?) and the coefficients in the action equ.(p|) as in the (j)^ case. 



4 Vertex Models on Planar (Fat) Graphs 

Writing down vertex model actions for the N 00 limit is a simple matter of putting a Tr in front of the 
previously discussed actions such as equ. (^ and interpreting the variables as N x N matrices rather than 
scalars. In the case of higher than 0"^ potentials the non-commutativity of the matrix variables will also 
allow the introduction of new terms compared with the thin graphs as, for example, terms such as [$^$]^ 
and [$'l']2$2^ where the $, are now NxN matrices, may be distinguished. This difference from the thin 
graph case arises because a consistent cyclic order of legs round a vertex may be defined in the planar limit 
which makes such distinctions meaningful, as can be seen in Fig. 4. At the level of the partition function 
an explicit contour integral to pick out the 27ith order is no longer necessary as the number of planar 
graphs increases only exponentially with the order, so the vertex coupling ( "cosmological constant") may 
be tuned to get a diverging number of vertices. This also means that the factorial prefactor Nn need no 
longer be peeled off. Writing down such partition functions is one thing, solving them is another as one 
is no longer dealing with saddle points in simple scalar integrals. In general the result is an insoluble, 
or at least unsolved, matrix model. Nonetheless, the planar graph vertex models can be mapped onto 
other, solved, models in some particular cases which we discuss below. 
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The non-commutativity of the matrix variables in the planar fat graph limit iV — > oo plays no role on 
4>^ graphs, as the cyclic symmetry of the trace ensures that possible potential terms are identical to the 
thin graph case. We have, writing ^ = X + lY for the N x N complex matrix $ and N x N Hermitian 
matrices X, Y 



Tr [($t)3 _,_ $3j ^ 2 Tr [X^ 
i Tr [$3 - i^^f] = 2 Tr [Y^ - 



-XY"^] 
YX^] 

- 3XY^] 
3YX^] 



(33) 



just as for thin graphs. The fat graph version of the action in equs.(|^,^) therefore differs only in the 
presence of a trace 



S 



Trr-{X^ 



Y' 



(« + 37) ^^3 



{8-^)X'^Y 



(34) 



Equ.(Q) is known to represent the Ising model on dynamical triangulations (i.e. the dual of planar 0^ 
graphs) with the spins situated at the triangle vertices, when /3 = (5 = and (7 — a)/(a + 87) = c* 
pC| ]. Its critical behaviour is therefore that of the Ising model coupled to 2D gravity, giving the KPZ 
critical exponents [|l2| rather than the flat 2D lattice Onsager exponents. We can say little about the 
more general case, as the solution of the resulting matrix model is not known. 

For fat (\y' graphs the non-commutativity does play a role as there are now more quartic invariants 
than in the thin case. Writing ^ = X ^ lY again, we have the following possibilities 



Tr [($t)2$2 

Tr ((4>^)$)^ 
Tr [($t)4 + $4j 
i Tr - ($t)4] 
Tr [($t)3$ + $3(^^t)] 
iTr [($t)3$_$3(^t)] 



Tr {X^ + + 2XYXY\ 

Tr {X^ + + 4X2y2 - 2XYXY\ 

2 Tr [X* + y* - 4X2y2 - 2XYXY\ 

8 Tr {XY^ - X^Y\ 

2 Tr {X^ - Y^\ 

4 Tr {X^Y + Y^X\ . 



(35) 



Dropping the terms containing an asymmetric combination of it is possible to construct an action 

that is equivalent to the Ising action for spins on the vertices of a random quadrangulation (i.e. the duals 
of graphs) by choosing the appropriate combination of the remaining terms 



Tr<\x^ 



iy2 _ 1^4 _ £!(xyxy + 2x2y2) - 



y^ 



(36) 



where c* is again the dual to to the Ising couphng c = cxp(— 2/3), c* — (1 — c)/(l + c)l That eq u.jjs^ ) 
dots represent an Ising model can be shown by retracing the steps that lead from equ.(|29|) to equ. (|30|) in 
the thin graph case, or by deriving the weights in equ.(|36|) from direct consideration of the Ising model 
on the random quadrangulation |10| . The most general 0'* vertex model containing a linear combination 
of the terms in equ.(|3^) remains unsolved, but some further examples can be solved, or at least reduced 
to eigenvalue integrals and solved in principle. 



One of these was noted some years ago in [13 by Ginsparg. If we consider the action 



S' = Tr|i$t$_ct 
12 

the resulting Hermitian matrix model action is 



($t$)%($t)2$2 



5* = Tr I i + y2) - 2a (X^ + y^ 



(37) 



(38) 



7 



which can be seen to be an 0{2)l XY model by introducing a further Hermitian matrix M and rewriting 
the potential term to give 

£• = Tr I i (X^ + y2 A/2) - 2^/^{X'^ + Y^)M^ . (39) 

Allowing oneself the luxury of two couplings rather than the one of the XY model also gives rise to 
an interesting class of models. The combination 



(40) 



which is equivalent to the real action 



5 = Tr I i (X^ + Y'^) - (a + /?) (X* + y^) - 2(a - I3)XYXY - AfSX^Y^ \ (41) 

can be shown to be a free fermion model |l4[ . 
Taking another combination 



gives 



5 = Tr |i (X^ + Y^) -{a + l3){X'^ + Y"^) - 2{a - f3)XYXY + AjSX^Y^' 

which is similar to the free fermion model apart from the sign of the X^Y^ coupling. 
If we take /3 negative and a = —j3 in equation. (|3|) we get the action 

S = Tri^{X^ + Y"^) - m3\{XYXY + X^Y^ 



(42) 



(43) 



(44) 

considered by Cicuta et.al. [ p^ in their discussion of a three-edge colouring problem for planar graphs. 
Although not yet solved exactly the partition function for this model can be reduced to an eigenvalue 
integral and its properties explored 



Z = jWdK WiK - \,f cxp (-\Y. An n (1 - (A, + A,)2) 



(45) 



where the Ai are the N eigenvalues of the N x N matrices. The scaling of the edge of the eigenvalue 
distribution of this model in the N oo limit suggest that the critical behaviour is Ising-like, as the 
exponents are identical. 

The same eigenvalue expression arises on taking a — f3 {(3 positive) in equation. (Bsh to give the action 



= Tr |i (X^ + Y^) - 2I3{X^ - Y^ f 



Adopting similar tactics to the XY model and introducing a further matrix M gives 



S = Tr\^{X^ + Y'^ + Af 2) - 2 V/3(X^ - Y^)^!^ 



(46) 



(47) 



or 

S = Tr I^Af^ + 1^2(1 - 2y^M) + ]^Y'^{1 + 2y^M)^ (48) 
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whose partition function reduces to the eigenvalue integral 

Z = JUdX. UiX, - \,f exp [ -i ^ An Y[{1 ~ 2^{\, + A,))-^(l + 2V^(A. + A,))-^ (49) 

i i<j \ i / ij 

which is identical to equation. (^) when the two measure factors are combined. This is not so surprising 
when one realizes that the actions of equ. (^) and equ.(|4^) are related by our old friend the linear 
transformation X ^ {X + Y) / V2, Y ^ {X - Y)/V2. 

We note in closing that the appropriate linear combination of the first three terms in equ.(|35|) produces 
an XYXY potential alone, which was also shown in to display Ising critical behaviour by carrying 
out an ingenious non-linear transformation on the matrices. Similar models occur in the enumeration of 
meanders Jl7[ . 



5 Discussion 



We have considered vertex models on thin (p^ and 0^ graphs and shown that various sorts of critical 
behaviour - Ising, Potts, multi-spin interaction, . . . may be obtained by restricting the most general 
couplings allowed by the arrow reversal symmetry in the vertex models to various subsets. In this sense 
the vertex models discussed here play a similar role to the eight vertex model on the regular square lattice, 
in whose phase diagram lurks the critical behaviour of many other familiar models One important 
difference from the regular lattice case is that the notion of orientational order for edges round a vertex is 
lost on random graphs (there is no "up" and "down"). We showed that writing the vertex model actions 
in real form allowed us to recognise transformations from other known models and deduce the critical 
behaviour. The simplicity of the thin graph models also allowed the explicit solution of the saddle point 
equations in some cases to confirm these deductions, in particular for the various Ising models. We noted 
that the vertex models, when written in their equivalent real form, also admitted a physical interpretation 
in terms of spin models on cacti graphs, or loop/line gases on the original random graphs. 

Although just as easy to write down actions for fat graphs, solving such models is a rather more 
difficult task than on thin graphs. We presented some solutions for vertex models on fat graphs by 
showing their equivalence to already solved models, mostly by following a similar tack to the thin graph 
exposition and rewriting the models in terms of two Hermitian matrices. In these examples Ising, XY and 
Ising-like (for the model of [p^) critical behaviour was observed, as well as a model that was equivalent 
to free fermions (that of equ.([4l|), jlj]). On fat graphs the cyclic order of edges round a vertex may be 
consistently defined, so the number of distinguishable terms in higher than <p'^ potentials is greater than 
in the thin graph case. 

As noted in the introduction, the idea of using integrals over complex variables and saddle point 
methods to extract the asymptotic behaviour of directed graph models for large graphs has been suggested 
by Whittle From a graph theory point of view all of the calculations in this paper are exercises in 
counting various sorts of more or less esoteric directed random graphs. The use of scalar "actions" and 
saddle point methods seems to be as effective in investigating the statistical mechanics of vertex models 
on thin graphs (or, if one prefers, the statistical mechanics of directed graphs) as similar methods have 
been for other statistical mechanical, notably spin, models living on thin graphs. 

DAJ would like to thank Thordur Jonsson for various discussions about vertex models on planar 
graphs. 
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Figure 2: A section of a Husimi tree composed of triangles. The underlying (j)^ graph is also shown. 
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Figure 3: A section of a Husimi tree composed of squares. The underlying cp"^ graph is also shown. 
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(a) (b) 



Figure 4: Two different cyclic orientations of a quartic vertex which may be distinguished in a planar 
graph: (a) [3>t$]2, (b) [$t]2$2 
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